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Statistical properties of the critical current distribution in superconductor with fractal clusters of 
a normal phase are considered. It is found that there is the range of fractal dimensions in which the 
variance and expectation for this distribution increases infinitely. Simple technique of avoiding such 
a divergence by the use of truncated distributions is proposed. It is suggested that the most current- 
carrying capability of a superconductor can be achieved by modifying the cluster area distribution 
in such a way that the regime of giant variance of critical currents will be realized. 



I. INTRODUCTION 

Considerable recent attention is being drawn to the fractal behavior of magnetic flux in type-II superconductors 
[1]- [3]. High temperature superconductors (HTS) containing clusters of correlated defects [4,5] are of special interest 
in this field. The case of clusters with fractal boundaries provides new possibilities for increasing the critical current 
value [2,6]. By virtue of the capability to trap a magnetic flux, such clusters can appreciably modify magnetic and 
transport properties of superconductors [3,7,8]. The distribution of the critical currents in superconductor containing 
the normal phase clusters with fractal boundaries has unusual statistical properties, and it is these features that will 
interest us. 

Let us consider a superconductor containing columnar inclusions of a normal phase, which are out of contact with 
one another. These inclusions may be formed by the fragments of different chemical compositions, as well as by 
the domains of the reduced superconducting order parameter. The similar columnar defects can readily be created 
during the film growth process [5,9,10]. In the course of the cooling below the critical temperature in the magnetic 
field ("field-cooling" regime) the magnetic flux will be trapped in the isolated clusters of a normal phase so the two- 
dimensional distribution of the flux will be created in such a superconducting structure. When the transport current 
is passed transversely to the magnetic field, this one is added to all the persistent currents, which circulate around 
the normal phase clusters and keep the trapped magnetic flux to be unchanged. By the cluster we mean a set of the 
columnar defects, which are united by the common trapped flux and are surrounded by the superconducting phase. 
Inasmuch as the distribution of the trapped magnetic flux is two-dimensional, instead of dealing with an extended 
object, which indeed the normal phase cluster is, we will consider its cross-section by the plane carrying a transport 
current. As was first found in Ref. [2], clusters of a normal phase can have fractal boundaries, and this feature has a 
significant effect on the dynamics of the trapped magnetic flux [3,8]. In the subsequent consideration we will suppose 
that the characteristic sizes of the normal phase clusters far exceed both the coherence length and the penetration 
depth. This assumption agrees well with the data on the cluster structure in YBCO films [2,3,10], as well as will allow 
us to highlight the role played by the cluster boundary in the magnetic flux trapping. 



II. GIANT DISPERSION OF CRITICAL CURRENTS 

Suppose that there is a superconducting percolation cluster in the plane of the film where a transport current flows. 
Such a structure provides for an effective pinning, because the magnetic flux is locked in finite clusters of a normal 
phase. When the transport current is increased, the trapped magnetic flux remains unchanged until the vortices start 
to break away from the clusters of pinning force weaker than the Lorentz force created by the current. As this takes 
place, the vortices must cross the surrounding superconducting space, and they will first do that through the weak 
links, which connect the normal phase clusters between themselves. Such weak links form easily in HTS characterized 
by an extremely short coherence length. Diverse structural defects, which would simply cause an additional scattering 
at long coherence length, give rise to the weak links in HTS. Weak links arise readily on twin boundaries, and magnetic 
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flux can easily move along them [11]. Whatever the microscopic nature of weak links could be, they form the channels 
for vortex transport. Accordingly to weak link configuration each normal phase cluster has its own value of the critical 
current, which contributes to the total distribution. By the critical current of the cluster we mean the current of 
depinning, that is to say, such a current at which the magnetic flux ceases to be held inside the cluster of a normal 
phase. The critical current distribution is related to the cluster area distribution, because the cluster of a larger size 
has more weak links over its boundary with the surrounding superconducting space, and thus the smaller current of 
depinning [3]. 

In the practically important case of YBCO films with columnar defects the exponential distribution of the cluster 
areas can be realized [2] , which is the special case of gamma distribution. The exponential distribution has only one 
characteristic parameter (mean cluster size) , so there are not many possibilities to modify the geometric morphological 
properties of the clusters in this simplest case. By contrast, in the case of gamma distribution there is an additional 
way for optimizing the cluster structure of the composite superconductor by the control of two independent parameters 
in the course of the film growth. One of the aims of the present work is to find how the cluster area distribution 
should be optimized in order to get the highest current-carrying capability of a superconductor. In Ref. [7] the critical 
current distribution was derived in the case of gamma distribution of fractal clusters areas, which has the following 
probability density 

A° exp(-A/A ) 

w (A) = xi - (1) 

r(. 9 + i)^ +1 

where T(v) is Euler gamma function, A is the cluster area, Ao > and g > —1 are the parameters of gamma 
distribution that control the mean area of the cluster A = (g + l)Ao and its variance a\ = (g + 1) Aq. For further 
consideration it is convenient to introduce the dimensionless area of the cluster a = A/ A, for which the distribution 
function of Eq. (1) can be rewritten as: 

w {a) = rfg+T) 1 ^ exp ( ~ {g + 1] a) (2) 

The mean dimensionless area of the cluster is equal to unity, whereas the variance is determined by one parameter only: 
<j\ = 1/ (g + 1). The probability density of Eq. (2) is presented in Fig. 1 for the characteristic values of g-parameter. 
The critical current distribution has the following form [7] 

m = r^JOH-D-i cxp (3) 

where 

G = 



69+ 1 - (£>/2) exp (0) + 1, 



D 

° = 9+l + ^ 

T(v,z) is the complementary incomplete gamma function, i = I/I c is the dimensionless electric current, I c = 
a (AqG) D ^ 2 is this the critical current of the resistive transition, a is the form factor, and D is the fractal dimension 
of the cluster perimeter. The value of D specifies the scaling relation P X I D oc A x l 2 between perimeter P and area A 
of the cluster [12,13]. 

The probability density curve for critical current distribution of Eq. (3) has the skew bell-shaped form with inherent 
broad "tail" extended over the region of high currents (see Fig. 2). As may be seen from this graph, the critical current 
distribution spreads out with some shift to the right as ^-parameter decreases. This broadening can be characterized 
by the standard deviation of critical currents 



_ » Tjg + l-D) / r( g + l-£>/2) V 

y r ( , + i) -{ r( 9 + i) ) (4) 

The standard deviation grows nonlinearly with increase in the fractal dimension, as is illustrated in Fig. 3. The 
peculiarity of the distribution of Eq. (3) is that its variance becomes infinite in the range of fractal dimensions 
D > g + 1. The distributions with divergent variance are known in probability theory - the classic example of that 
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kind is Cauchy distribution [14]. However, such an anomalous feature of exponential- hyperbolic distribution of Eq. (3) 
is of special interest, inasmuch as the current-carrying capability of a superconductor would be expected to increase in 
the region of giant variance. Then the statistical distribution of critical currents has a very elongated "tail" containing 
the contributions from the clusters of the highest depinning currents. 

The distribution of Eq. (3) has one even more striking feature: its mathematical expectation, which represents the 
mean critical current 

is also divergent in the range of fractal dimensions D ^ 2 {g + 1). At the same time, the mode of the distribution 

modef(i) = (G/9) D ^ 2 remains finite for all possible values of the fractal dimension 1 ^ D ^ 2. 

The reason for divergence of the mean critical current consists in the behavior of the cluster area distribution of 
Eq. (2). As may be seen from Fig. 1, the graph of the distribution function of Eq. (2) takes essentially different shapes 
depending on the sign of (^-parameter - from the skew unimodal curve (1) (at g > 0) to the monotonic curve (3) with 
hyperbolic singularity at zero point (at g < 0). In the borderline case of g — 0, which separates these different kinds 
of the functions, the distribution obeys an exponential law (curve (2)). It is just for negative values of ^-parameter 
that the mean critical current diverges. The contribution from the clusters of small area to the overall distribution 
grows at g < 0. Since the clusters of small size have the least number of weak links over a perimeter, they can best 
trap the magnetic flux. Therefore, an increase of the part of small clusters in the area distribution of Eq. (2) leads to 
a growth of the contribution with high depinning currents made by these clusters in the critical current distribution 
of Eq. (3). Just as a result of this feature the mean critical current diverges at g < 0. Nevertheless, the total area 
between the curve of the probability density (like curve (3) in Fig. 2) and the abscissa remains finite by virtue of the 
normalization requirement. 



III. TRUNCATED DISTRIBUTION OF CRITICAL CURRENTS 



Obviously, the proper critical current cannot be infinitely high as well as the clusters of infinitesimal area do not 
really exist. There is the minimum value of the normal phase cluster area A m , which is limited by the processes of 
the film growth. So in YBCO based composites, which were prepared by magnetron sputtering on sapphire substrates 
with a cerium oxide buffer sublayer [3] , the sample value of minimum area of the normal phase cluster has been equal 
to A rn = 2070 nm 2 at mean cluster area A = 76500 nm 2 , that corresponds to the lower bound of the dimensionlcss 
area of the cluster a m = A m j A — 0.027. In view of this limitation, we will describe the distribution of the normal 
phase cluster areas by the truncated version of the probability density of Eq. (2): 

w(a\a^ dm) = — zw(a) (6) 

1-W (dm) 

where 7 (v, z) is the complementary gamma function, h (x) = | q 'j,^ ^ < q * S ^ e Heaviside step function, 

Am 

W (a m ) = J w(a)da= T (g + 1) ^ ' 



is the truncation degree, which is equal to the probability Pr{Va < a m } to find the cluster of area smaller than the 
least possible value a m in the initial population. 

The expression of Eq. (6) gives the conditional distribution of probability, for which all the events of finding a 
cluster of area less than a m are excluded. Thus the truncation provides a natural way to fulfil our initial assumption 
that the cluster size has to be greater than the coherence and penetration lengths. New distribution of cluster areas 
gives rise to the truncated distribution of the critical currents: 

where i m = (G/ (g + 1) a m ) D ' 2 is the upper bound of the depinning current, which corresponds to the cluster of the 
least possible 

Then, instead of Eqs. (4) and (5), the standard deviation and mean critical current are 
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I = 



r(g + l-D/2,(g + l)a m ) c n 
T(g + l,(g + l)a m ) 



For the truncated distribution of Eq. (8) the possible values of depinning currents are bounded from above by 
the quantity i m , therefore the mean critical current as well as the variance do not diverge any more. Both of 
these characteristics are finite for any fractal dimensions, including the case of maximum fractality (D = 2). The 
corresponding graphs for standard deviation are presented in Fig. 4. All the curves are calculated at g = —0.2 (as 
for the main curve (2) of the mean critical current in Fig. 5). In this case the variance for initial distribution of the 
critical currents is infinite, so no graph for W (a m ) = is shown in Fig. 4 at all. The dependence of the mean critical 
current on the fractal dimension in the case of truncated distribution is demonstrated in Fig. 6. The corresponding 
graphs in Figs. 4 and 6 are drawn for the same values of truncation degree. The truncation degree is related to the 
least possible area of the cluster by the equation (7). The values of W (a m ) and a m involved in Figs. 4 and 6 are 
presented in the Table I. 

The degree of truncation gives the probability measure of the number of normal phase clusters that have the area 
smaller than the lower bound a m in the initial distribution. If the magnitude of W (a m ) is sufficiently small (no more 
than several percent), then the procedure of truncation scarcely affects the initial shape of the distribution and still 
enables the contribution from the clusters of zero area (therefore, of infinitely high depinning current) to be eliminated. 
It is interesting to note that the very similar situation occurs in analyzing the statistics of the areas of fractal islands 
in the ocean [12]. The island areas obey the Pareto distribution, which also has the hyperbolic singularity at zero 
point that causes certain of its moments to diverge. For exponential distribution of the cluster areas, which is valid 
in the above-mentioned case of YBCO films [3], the probability to find the cluster of area smaller than a m in the 
sampling is equal to Pr{Va < a m } = 2.7% only. In principle, the truncation procedure could be made here, too, but 
there is no need for that, because the contribution of infinitesimally small clusters is finite at g — (and equal to zero 
at g > 0). 

Referring to Figs. 4 and 6 it can be seen that the dependencies of the standard deviation and the mean critical 
current on the fractal dimension become smoother and smoother with increase in the degree of truncation. At the 
same time, the inherent tendency of initial distribution is still retained: as the fractal dimension increases, the critical 
current distribution broadens out (the variance grows, see Fig. 4), moving towards higher currents (the mean critical 
current grows, see Fig. 6). As may be seen from Fig. 2, this trend is further enhanced with decreasing g-parameter. 
The most current-carrying capability of a superconductor should be achieved when the clusters of small size, which 
have the highest currents of depinning, contribute maximally to the overall distribution of the critical currents. Such a 
situation takes place just in the region of giant variance of critical currents. So far, the least magnitude of (/-parameter 
(equal to zero) has been realized in YBCO composites containing normal phase clusters of fractal dimension D = 1.44 
[2]. The critical currents of superconducting films with such clusters are higher than usual [6,9,10]. It would thus 
be expected to further improve the current-carrying capability in superconductors containing normal phase clusters, 
which will be characterized by area distribution with negative magnitudes of ^-parameter, especially at the most 
values of fractal dimensions. 

Thus, it has been revealed that the distribution of the critical currents in superconductor with fractal clusters of 
a normal phase has the anomalous statistical properties, implying that its variance and expectation diverges in the 
certain range of fractal dimensions. It may be expected that the most current-carrying capability of a superconductor 
can be achieved by optimization of the cluster area distribution that involves reducing g-parameter with concurrent 
increasing the fractal dimension. 
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TABLE I. Degree of truncation of the cluster area distribution and corresponding least possible area of the normal phase 
cluster at g = - 0.2 



Truncation degree W (a m ) 


Minimum area a m 


Markers in Figs. 4, 6 


0.01% 


1.144 x 10"° 


a 


0.1% 


2.034 x 10~ 4 


b 


1% 


3.622 x 10~ 3 


c 
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FIG. 1. The distribution of the areas of normal phase clusters at different values of g-parameter. Curve (1) corresponds to 
the case of g = 1; curve (2) of g = 0; curve (3) of g — —0.5. The arrow indicates the mean cluster area. 
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FIG. 3. Influence of the fractal dimension of the cluster boundary on the standard deviation of critical currents. Curve (1) 
corresponds to the case of g = 0.25; curve (2) of g = 0.5; curve (3) of g = 0.75; curve (4) of g = 1. 
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FIG. 4. Standard deviation graphs for truncated distribution of the critical currents at g = —0.2 with the different degree of 
truncation: curve (a) is drawn for W (a m ) = 0.01%; curve (b) for W (a m ) — 0.1%; curve (c) for W (a m ) = 1%. 



9 




1.0 1.2 1.4 1.6 1.8 2.0 



D 



FIG. 5. Influence of the fractal dimension of the cluster boundary on the mean critical current. Curve (1) corresponds 
the case of g = —0.4; curve (2) of g = —0.2; curve (3) of g = 0; curve (4) of g — 1. 
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FIG. 6. Mean critical current graphs for truncated distribution at g = —0.2 with the different degree of truncation: curve 
(2) is drawn for W (o m ) = 0; curve (a) for W (o m ) = 0.01%; curve (b) for W (a m ) = 0.1%; curve (c) for W (a m ) = 1%. The 
dotted lines are given for truncated distribution, the solid line for initial distribution. 
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